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Using a time-dependent linear (Rabi) coupling between the components of a weakly interacting 
multicomponent Bose-Einstein condensate (BEC), we propose a protocol for transferring the wave- 
function of one component to the other. This "Rabi switch" can be generated in a binary BEC 
mixture by an electromagnetic field between the two components, typically two hyperfine states. 
When the wavefunction to be transfered is - at a given time - a stationary state of the multicom- 
ponent Hamiltonian, then, after a time delay (depending on the Rabi frequency), it is possible to 
have the same wavefunction on the other condensate. The Rabi switch can be used to transfer 
also moving bright matter-wave solitons, as well as vortices and vortex lattices in two-dimensional 
condensates. The efficiency of the proposed switch is shown to be 100 % when inter-species and 
intra-species interaction strengths are equal. The deviations from equal interaction strengths are 
analyzed within a two-mode model and the dependence of the efficiency on the interaction strengths 
and on the presence of external potentials is examined in both ID and 2D settings. 



I. INTRODUCTION 



The past decade has witnessed a tremendous explosion of interest in the experimental and theoretical studies of 
Bose-Einstein condensates (BECs) P, [2| • Numerous aspects of this novel and experimentally accessible form of matter 
have been since then intensely studied; one of them concerns the investigation of the behavior of multicomponent 
BECs, which have been experimentally studied in either mixtures of different spin states of ^'^Na [1, 0, Q or *^Rb 
UM, i, M, [U [H, [iMl M, M,M,M,M,^, or even in mixtures of different atomic species such as '^^K-^'^Rh 
and ^Li-i33Cs [237^ 

The dynamics of a multicomponent BEC is described, at the mean-field level, by coupled Gross-Pitaevskii (CP) 
equations, taking into account the self- and cross- interactions between the species. In this framework, a number of 
properties and interesting phenomena have already been extensively analyzed. Among them, one can list ground state 
solutions [HjH^,!!^ and small-amplitude excitations [13] of the order parameters in multicomponent BECs, as well as 
the formation of domain walls [1^ and various types of matter- wave soliton complexes [1^ , spatially periodic states [s^l 
and modulated amplitude waves [sij . Quantum phase transitions in Bose-Bose mixtures have been investigated both 
theoretically [H, [H, [H, [3^ and experimentally ^221]. Moreover, several relevant works analyzed different aspects of 
purely spinor (F = 1) condensates (which have been created in the experiments [ 3 | . [l3 l ), i ncluding the formation of spin 
textures [3 , sp in domains [s^, various types of vector matter- wave solitons [371 l38l, [sol l40j . studies of ferromagnetic 
properties [41|, and so on. 

An important resource for the experimental control of multicomponent BECs is the possibility to use a two-photon 
transition to transfer an arbitrary fraction of atoms from one component to another, e.g., from the |1, — 1) spin state 
of *^Rb to the |2, 1) state. The transfer can also occur by using an electromagnetic field inducing a linear coupling, 
proportional to the Rabi frequency, between the different components. In Ref. [soj it was shown that, in analogy 
with systems arising in the field of nonlinear fiber optics (such as a twisted fiber with two linear polarizations, or an 
elliptically deformed fiber with circular polarizations [i^l), exact Rabi oscillations between two condensates can be 
analytically found when inter-species coupling are equal to unity (in proper dimensionless units). 

In this paper, we propose a protocol enabling the transfer of the wavefunction of a condensate to another, even 
in presence of interactions. The proposed protocol requires a time-dependent Rabi frequency: this "Rabi switch" is 
realized by turning-on the linear coupling for a pertinent period of time, so as to transfer the maximal fraction of the 
condensate from the first to the second component. The efficiency of the switch is maximal, if all interaction strengths 
(nonlinearity coefficients) are equal. If one deviates from the ideal case, the efficiency is modified: to analyze more 
realistic situations, we show that it is possible to effectively describe the deviations from equal interaction strengths by 
a two-mode ansatz, where the impossibility of transferring all the particles from a condensate to the other is identified 
as the self-trapping of the initial condensate wavefunction. Even though in the original experiments (see e.g. [^) the 
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Rabi coupling was used to transfer ground states between two repulsive condensates, our protocol can be efRciently 
used for transferring also matter-wave solitons in one-dimensional (ID) attractive condensates, as well as vortices and 
even vortex lattices in two-dimensional (2D) repulsive condensates. We study the efficiency of the proposed Rabi 
switch in each of these situations and we discuss the generalization of the same idea to a 3-component condensate, 
where our approach would realize a "Rabi router" of matter into desired components. 

The protocol proposed in this paper would allow for to copy a wavefunction from a condensate to the other in 
the presence of either attractive or repulsive interactions among atoms, and this could improve the efficiency in 
the experimental manipulation of matter solitons and vortices; from this point of view, it provides a matter-wave 
counterpart for optical switches realized in nonlinear fiber optics, which are important tools to control optical solitons 

0- . . 

Our presentation is structured as follows. In Section II we present the theoretical framework needed to describe 
the Rabi switch in a two-component Bose gas, which is valid for attractive or repulsive interactions. In Section III 
the generalization to multicomponent BECs is discussed. In Sections IV and V we provide the results of our analysis 
in ID and 2D settings, respectively; there, we also show how the external trapping potentials affect the efficiency of 
the Rabi switch, and compare the findings of the two- mode model with numerical results. Finally, in Section VI we 
present the conclusions and outlook of this work. 



II. THE RABI SWITCH 



The prototypical system we consider is a two-component Bose gas in an external trapping potential: typically the 
condensates are different Zeeman levels of alkali atoms like ^^Rb. Experiments with a two-component ^^Rb condensate 
use atom states customarily denoted by |1) and |2); in particular, the states can be \F = 2,mp = 1) and |2, 2), like, 
e.g., in or 1) and |2,1), like, e.g., in [1] (see also the recent work [1^). In general, the condensates |1) and 
1 2) have different magnetic moments: then in a magnetic trap they can be subjected to different magnetic potentials 
Vi and V2, eventually centered at different positions and having the same frequencies (like in the setup described in 
Q) or different frequencies [l^. In the ratio of the freq uencies of V2 and Vi is v^. It is also possible to add a 
periodic potential acting on the two-component Bose gas [ll|, [22| . 

The two Zeeman states |1) and |2) can be coupled by an electromagnetic field with frequency tJoxt and strength 
characterized by the Rabi frequency fi/j, as schematically shown in Fig.^ A discussion of (and references on) the 
experimental manipulation of multicomponent Bose gases can be found in [4J,|45| (for a recent experimental realization 
of this coupling, see also [20l|). The detuning is defined as Wcxt — '-^Oi where fvjjQ is the energy splitting between the two 
states (e.g., in [1C| ~ 27r x 2 MHz). For concreteness we assume that the Rabi coupling can be turned on starting 
at a given instant, say > 0; at later times, the Rabi coupling coherently transfers particles between |1) and |2) at a 
Rabi frequency Qn. When the number of components is larger than two, more coupling electromagnetic fields could 
similarly be added. The transfer of particles between hyperfine levels may be also regarded as an "internal Josephson 
effect" ■ since it is similar to the Josephson tunneling of particles between Bose condensates in a double-well potential 
[46l |47| ; the only difference is that in the "internal Josephson effect" the two condensates are spatially overlapping, 
while the left and right part of a single-species BEG in a double-well are separated by the energy barrier. Thus, the 
roles of the Rabi frequency and the detuning in the internal Josephson effect are analogous to the ones played by the 
tunneling rate and the difference between zero-point energies of the two wells, respectively. 

In the rotating wave approximation, the dynamics of the two-component Bose-Einstein condensates is described by 
two coupled GP equations [1, |H, E^, which, in a general 3D setup and in dimensionless units, read 



.d'4'2 



1 



A-l-yi(f) -l-giilViP +ffi2|V'2| 



1 



A + V2(r)+5l2|V'l|^+522|V'2| 



■01 + a{t)'4>2, 
■02 + a{t)ipi, 



(1) 

(2) 



where ^j(f,t) are the wavefunctions of the j-th condensate {j — 1,2), Vj are the respective trapping potentials 
(typically, an harmonic potential and/or an optical lattice, plus the eventual detuning, absorbed in them) and the 
quantities gij , which are proportional to the scattering lengths a.^ of the interactions between the species i and the 
species j, describe the intra- {j — i) and inter- (j ^ i) species interactions. The system ([Ij-© consists of two linearly 
and nonlinearly coupled GP equations: the linear coupling is provided by the Rabi field, while the nonlinear coupling 
is proportional to gi2 and is due to the scattering between particles of the different species. The scattering lengths 
for experiments with Zeeman levels of ^^Rb atoms are normally quite similar: in fact, if |1) = |2, 1) and |2) ~ |1, —1) 
the scattering length ratios are an : 012 : 022 = 0.97 : 1.00 : 1.03, while if |1) = |2, 1) and |2) = |2,2) they are 
an : ai2 : 022 = 1.00 : 1.00 : 0.97. Furthermore, one of the a^'s can be varied through Feshbach resonances [l|. 
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FIG. 1: Josephson coupling of the two Zeeman levels |1) and |2) through an electromagnetic field with frequency ix)oxt and 
strength characterized by the Rabi frequency Qr. 



The term a is proportional to the Rabi frequency ^Ib [50], and serves the purpose of transferring the condensate 
wavefunction of |1) to condensate |2) and of controUing the time modulation of the Rabi frequency. 

When the external potentials are the same (Vi — V2 — V) and the interactions strengths are equal {gn — 1722 — 
9i2 = 5)1 Eqs. (HI-© can be rewritten in a more compact form as: 

= -^At/; + i^^Gi^)^ + F(r)V + a{t)P^, (3) 

where 

As a result of the fact that G and P commute, one can decompose the solution ip of the "inhomogeneous" problem 
described by Eq. ^ as 

i,{r,t) = Um{r,t), (5) 
where U(t) is the matrix of the homogeneous problem 

um^^HPm^i^X,-:^^). (0) 

with X{t) ~ J* a{t')dt' . Substituting Eq. ([6]) in Eq. ([3]), it is readily found that 0(r, t) satisfies the evolution equation 

M 

i^ = ncP^-^Aq^+{cl,^Gc^)(t, + V{r)c^, (7) 

which is identical to Eq. ([3]), but without the Rabi term proportional to a{t)P. 

When the external potentials or the interaction strengths are different, it is formally possible, as discussed in the 
Appendix A, to perform a decomposition like the one given in Eq. ([5]) and remove the Rabi term. However, this is 
done at the price of introducing time- and space- dependent effective interaction strengths in the nonlinear terms: in 
particular, for different external potentials the effective interaction strengths are both time- and space- dependent, 
while when the interaction strengths are different, a nonlinear Josephson term also arises in the nonlinear terms (i.e., 
is proportional to (j)2 through terms proportional to products As a result, the removal of the Rabi term 

is of little practical use and one has to resort to numerical or variational estimates: in the following, we will focus 
on the situation of different g^ 's, and show that the efficiency of the Rabi switch for small deviations from the equal 
strengths situation can be effectively described by a two-mode model. 

When the Rabi frequency is fixed, oscillations of atoms between the two components have been studied [1, [l^, [s^, 
liol . [sil . [53 , [53I and experimentally observed [1, . In the following, we consider a time-dependent Rabi frequency 
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a{t): through a proper choice of a a{t), we can transfer the wavefunction of |1) to |2). More precisely, we propose 
a way to perform the following operation: at a time one has all the particles in |1) in the wavefunction ipi^VjiQ), 
and no particles in |2) {ip2{r,to) = 0). At a time ti we wish to have all the particles in |2) in the same wavefunction: 
'>p2i'f^, ti) = ipiir, to) (apart a phase factor). The protocol proposed in this paper allows the transfer of the wavefunction 
if ipi{r, to) is a stationary state (or a moving soliton, as discussed in Section IV) of the nonlinear Hamiltonian Ti. defined 
in ([7]). If it is not, we can however have at the time ti all the particles in |2) in the wavefmiction the condensate |1) 
would have had without the Rabi coupling [see Eq. ([13])] . 

The proposed protocol works with or without nonlinearity: but with the nonlinearity on, one can transfer also 
a soliton wavefunction; for instance, in ID one can have a matter- wave soliton of the species |1) propagating with 
velocity v, and after a time delay, the proposed Rabi switch will generate the same soliton with the same velocity 
in the species |2). Two remarks are due: (i) the transfer mechanism has the highest possible efficiency for equal 
interaction strengths, however it is very good and its efficiency is close to 1 in a wider region in the relevant parameter 
space [in ID, deviations from the integrable case of equal interaction strengths are discussed through a two-mode 
ansatz] ; (ii) the proposed mechanism is not copying the full many-body wavefunction of the weakly interacting Bose 
gas, but only the order parameter, which is related to the one-body density matrix. 

To be more specific, we assume that a{t) depends on time as 

< i < to, 

a{t) = <( 7, to<t<t,^to+6, (8) 

t > ti, 

where to is the switch-on time and 6 denotes the duration of the Rabi pulse. From Eq. ([5]) we readily find 

< t < <o, 

^it) = { lit -to), to<t<t,, (9) 

t > ti. 

Introducing the vector field by the decomposition in Eq. ([5]), it is observed that at the time to-, i.e., before the 
switch-on of the Rabi pulse, ^{r,to) = 4>{f,to), while for to <t <ti we find that 

-01 (r, t) = cos [-yit - to)]'/>i (r, t) - i sin [-/{t - to)](t)2{r, t), 
ip2{r,t) = -ism[-f{t-to)](l)iir,t) +cos[-f{t~to)]<p2ir,t). 



When the pulse duration is 



'5=f, (11) 
27 



at the end of the pulse we obtain 

'^i{r,ti) = ~i(j)2{r,ti), 



Mr,ti)^~tMr,ti). ^^^^ 

Since, in the interval [tojii], satisfies the homogeneous coupled GP equations ([7]), i.e. the same equations satisfied 
by the vector field ■0 in the interval [0,to], then 

V'2(^tl) \ ^ _ . iH(ti-to) f Mr, to) 



i^i{r,ti) J \ ip2{r,to) J 

Notice that if, instead of Eq. ([5]), one allows for a different time dependence of a{t), the time ti at which Eq. (fT^]) 
holds is given by the condition cosl(ti) — 0. For instance, with a{t) = for t < to and t > to + S, and a{t) = f{t) 
for to < t < to + 5, the pulse duration S such that Eq. (fT5|) is valid is given by dt' f(t' — to) = Tr/2. 

We are interested in the situation in which no particles are in |2) at to ('!/'2(^, to)=0); this situation may occur, e.g., 
in the case where only a single condensate has been prepared (if, eventually, particles exist in the other component, 
it is possible to remove them by the suitable application of a Rabi pulse). Notice that this has been experimentally 
realized e.g. in the experiments of (see also references therein). In such a case, Eq. implies that at the 

end of the pulse one has that (apart from a phase factor) the wavefunction describing the condensate |2) is the same 
wavefunction which the condensate |1) would have had in ti in the absence of the Rabi pulse. We can quantify the 
success of the described protocol in several ways: one of them will be to define the "efficiency" T as the fraction of 
atoms we are able to transfer from |1) to |2), i.e., 

^2(^1) . . 
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where Ni{t) = J df\ipi {r,t) p is the number of particles in the condensate i at time t. Notice that such a definition 
can even be extended in cases where the number of atoms in the second component is not zero initiaUy by replacing 
in the numerator of Eq. (|14p N2{ti) — > (A^2(^i) — ^2(^0))- Another more stringent way is to define a kind of "fidelity" 
F of the wavefunction transfer, i.e., the quantity 

J dr\r2ir,h)\ ■ \Mr,to)l (15) 

From Eq. (fT3)) we see that the efficiency of Rabi switch is 1 for equal interaction strengths, but the fidelity is not. 
However, we can have the fidelity to be equal to 1 if is a stationary state of Ti, corresponding to the eigenvalue /z, 
i.e. 

HtP^ir^ ^ ^i^^{r^. (16) 
If at time t = 0, i}{f, 0) = then 0(f, to) = e~^''*«-0^(f) and 

V'i(r,ii) = e"*'"' [cos(7(5) V'i(r,to) - « sin (7^) ?/'2(^ to)] , 



Mr,h) = e-'t'^ [-ism{-f5)Mr,to) + cosi-f5)Mr,to)] ■ ^^^^ 
When no particles are in |2) at to and the pulse duration is given by Eq. pT|) . one has 



(18) 



Equations ()18|) show that the wavefunctions of the two components have been exchanged up to a phase factor. This 
remarkable feature allows us, again with equal interaction strengths , to transfer the condensate wavefunction (with 
100% efficiency) from a populated hyperfine state to an empty one. In the following section we discuss how to transfer 
from a condensate to the other the wavefunction of a moving matter-wave soliton. 

We should further note about the latter that the nature of linear operators in the right hand-side of Eq. ([3]) is 
irrelevant in the derivation of Eq. ([5]). Hence, our analysis can be used to deal with: 

• repulsively interacting as well as attractively interacting systems; 

• continuum, as well as discrete systems; 

• homogeneous systems (in the absence of external potentials) or inhomogeneous systems (e.g., in the presence of 
external harmonic trap and/or optical lattice potentials); 

• one-dimensional systems or higher-dimensional ones. 

In what follows, we illustrate the versatility of the Rabi switch by examining characteristic examples for each of 
the above settings. We will then illustrate, how the perfect efficiency of the matter wave transfer (discussed above 
for equal inter-particle interactions) is "degraded" in more realistic situations (where such interactions are no longer 
equal) . 



III. GENERALIZATION TO MULTICOMPONENT BOSE-EINSTEIN CONDENSATES 

The protocol discussed in the previous Section can be generalized for Af > 2 components with a suitable choice of 
the time dependence of the Rabi frequencies atj transferring particles from the condensate i to the condensate j. For 
instance, for A/" = 3 and for equal potentials (Vi = V2 = V3 = V) and interaction strengths {gij = g), the relevant 
system of the three coupled GP equations can be written in the form of Eq. ([3]) , namely 

= -^Aij + {iP^GtP)tP + V{r)tP + Pit)tP (19) 

with 



/O ai3(t)\ 
V'= V'2 ,G = <? 010 ,P- ai2(t) a23(t) 

V V'3 / V 1 y V "23(t) J 



(20) 
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For general , the decomposition = U<j) with 



U = 



(21) 



fails to recast Eq. p9p in the homogeneous form 54]; however, it still removes the Rabi term when aij{t) = a{t) for 
any With iP = U<j> and U given by Eq. one finds i|f = -^A(f> +{((>'' G(t))(j). The matrix U,j{t) = 1,2,3) 

has diagonal elements Ujj = (1/3) [2 exp (iX) + exp {—2iT)] and off-diagonal ones Uij = Ujj — exp (iX). Once the Rabi 
term has been removed, the "Rabi switch" described in the previous Section can be applied also for general M to 
transfer a wavefunction from a condensate to any one of the others. 



Another choice of aij allowing for the removal of the Rabi term is provided by the generalization of Eq. 



namely 



0, 
0, 



0<t<to, 
to<t< ii. 



(22) 



with all the turned on/off at the same time, but with eventually different intensities. As an example, for J\f = 3, 
one may consider 712 = ai, 713 = a2 and 723 = 0: the matrix U(t), for to < t < ti, then reads 



/ riil + r2) airi 



r2) 

1— r2 



1 — r-2 



a2ri 



\ n2ri . " 



a2+airi(l+r2) 



aia2ri 



1— r2 



a^+a2 



aia2ri' 



a^+a2ri(l+r2) 



(23) 



where 



2ri 



.-1/2 



exp 



(24) 



This allows us to determine the transfer of matter from the first to the second and third component. Similar results 
may be obtained in the more general case of A/" > 2 components, i.e., "desirable" amounts of matter can be controUably 
directed to different hyperfine states according to their Rabi couplings. This general "Rabi router" is quite interesting 
in its own right, and it could be experimentally implemented in = 1 spinor condensates (sl. [Tsj . 



IV. RESULTS FOR ID SETTINGS 



We now consider the ID version of Eqs. ([T|)-(I21), which is relevant to the analysis of "cigar-shaped" condensates 
confined in highly anisotropic traps [H, [3| , and in dimensionless units reads 



.dip2 



+ Vi{x) -f .911 1^1 P + 9l2\4'2\ 



V2{x) + .912 I ■01 



.922|0'2p 



-01 -t- a{t)i>2, 
i^2 + a(OV'i- 



(25) 
(26) 



We use wavefunctions ■ipi{x,t) normalized to unity, so that Ni{t) + N2{t) = 1, with Ni{t) — J dx\tpi{x,t)\'^ . When 
the effect of the external potentials Vi is negligible (as, e.g., in the case of potentials varying slowly on the soliton 
scale) and in absence of the Rabi coupling (a = 0), the system p5|) -p6 )) becomes the Manakov system [5^, which 
is integrable for gn — gi2 — §22- In what follows we examine both attractive and repulsive interatomic interactions, 
corresponding, respectively, to negative and positive values of the scattering lengths, and we will consider the effect 
of the presence of the trapping potential on the wavefunction transfer. 



A. Stationary bright matter-wave solitons 



We consider in this subsection attractive interactions, gij < 0, in absence of external potentials. Putting iij = —gij, 
we first consider the ideal case where in = l\2 = I22 = ^ and assume that, at i = 0, all the particles of |1) are 
described by the 1-soliton solution of the nonlinear Schrodinger equation; thus, ip2ix, 0) = and 



^i(x,0) 



Vi/2 



cosh(£a;/2)' 



(27) 
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FIG. 2: Transferring a stationary bright matter-wave soliton: in (a)-(b)-(c) the density pj = is plotted for both components 
(|1) solid line; |2) dashed line) at the times t = to, to + 0.45, ti. In (d)-(e)-(f) the density is plotted at the same times, but for 
a velocity i; = 1 {S = n/2, ti = 5). 



In these units the soliton's chemical potential is ^ — — £^/8; turning on the Rabi coupling 7 at time Iq, and then 
turning it off at ti ^ to + 6, one gets, for to < t < ti, 

Ve/2 ^_,^(t-to) f cos7(i-io) 



' ' cosh(£a;/2) \ -ism-f{t - to). J ^ ' 

For 6 = 7r/(27) no particles are in the condensate |1) at ti, and the soliton wavefunction has been transferred in |2), 
i.e., ip2{x,ti) — —{i/2)e~^^^VI/ cos\i{ix/2). The transfer of the soliton wavefunction is illustrated in Fig. [2la)-(c). 
Notice that if we choose as initial condition 



(0) 



Mx,0) J cosh(V2) V V^2(0)e*'^2W 



(29) 



i.e., two bright solitons with particle numbers iVi(O), iV2(0) and phase difference A(/7(0) — (p2{0) ~ </'i(0), then, at 
t = ii, we obtain 

Ni{ti) = (cos (7(5) VA^i(0) + sin (-/S) sin Av3(0) ^^^2(0)) ^ + N2iO) sin^ {-fS) cos^ Aip{0). (30) 

This shows that, by choosing properly the pulse duration and the initial phase difference, one can transfer a "desired" 
part of the soliton wavefunction from one condensate to the other. 

Let us discuss now the interesting situation of different interaction strengthsithe aim there is to study the efficiency 
of the Rabi switch of the soliton wavefunction and qualitatively understand the effect of the deviation from the ideal 
case. To that effect, we introduce a variational two-mode ansatz and confine ourselves to the situation in which no 
particles are initially in |2). For to < t < ti, we choose the variational vectorial wavefunction 



_ , 4'vi{x,t) \ _ ( yiVi(i)e*'^^(*)<i>i(x) , 

- ' i^.2{x,t) [ y]V2(i)e^^^^*^^2(x) ' ' ^-^^^ 



where 



"^^"^^ = cos?fS/2) - (32) 

The variational parameters are the numbers of particles Ni(t) and their phases ipi(t). The variational vector wave- 
function (j3ip has been used in [56| to study the wavepacket dynamics for two linearly coupled nonlinear Schrodinger 

equations with £11 = £22 and ^12 = 0. For general iij's, the Lagrangian C = ~ ~5^'^v) — {^yTi-ipv) [where 

Ti. is given by Eq. (|B1[) and () denotes spatial integration], is computed in Appendix B, where we show that the 
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variational equations of motion for A''! — N2 and (^2 — fi are the equations of a (non-rigid) pendulum. The mass 
M of the pendulum depends on the ^^^'s according Eq. (|B7p and for In = £12 = £22 the pendulum mass is zero, 
allowing for the transfer of all the particles from a species to the other. When £11 ^ £22-, a detuning term in the 
pendulum equations is present [see Eq. (jBSp j. In the following, we will focus for simplicity on the more illmninating 
case £\\— £22-, with a general £\2- 
Introducing the variables 

r/ = iVi-iV2; <y3 = '^2 - <y5l, (33) 



one gets the equations of motion 



7) = 27^1 — 77^ sin 

Cp - -27^=2= cos + fll^i^T], 



(34) 



with initial conditions 77(^0) = 1 (i-e-j all the particles initially in |1)) and (/^(to) = 0. It is worth noticing that Eqs. 
([M)) are the same equations governing the tunneling of weakly-coupled BECs in a double- well potential i^Sj iSZ] (the 
only difference being that 7 corresponds to —K, where if > is the tunneling rate, which gives the same results 
for — > -|- 7r). Equations p4p are formally identical to the equations for an electron in a polarizable medium, 
where a polaron is formed f58| . Analytical solutions have been found for the discrete nonlinear Schrodinger equations 
describing the motion of the polaron between two sites of a dimer [H, [5^ . 

Eqs. ([34|) are the equations of a non-rigid pendulum 0, [s^l, with the effective Hamiltonian being 

Heff = Y?7^ + 27Vl-??2 cos^, (35) 



where the pendulum mass is given by 



M^^i /^' (36) 



When £11 — £12, then the mass in Eq. ([M]) vanishes and ij — —47^77. The duration S of the pulse needed to have a 
perfect switch is such that ri{ti ^ tQ + S) = —1, i.e. 6 = tt/ (27) in agreement with Eq. pTjl . If the mass M is positive 
(i.e., 512 < gii), then it is still possible to transfer all the particles from |1) to |2), provided that the mass M is smaller 
than a critical value Mc. If M < Mc, then the time ti at which ry(ti) — —1 will be different from -|- tt/{2j) (the 
analytical computation of the tunneling period for a mass M 7^ is reported in the Appendix of [57|). This means 
that for deviations from the ideal case, one can optimize the transfer by choosing a pulse duration different from Eq. 
(fTTjl . This is illustrated in Fig. [31 where we compare ri{t) obtained from the two-mode equations ([M]) with the results 
of the numerical solution of the GP equations ((25l) - ((26)) for £i2/-^ii = 13: the numerical and variational results are in 
good agreement for a wide range of the parameters (see the inset of Fig. [3]). The computation of Mc is done according 
to the method discussed in [4^ : namely, one has to compute the point at which self-trapping occurs, and the result is 

Me = 47. (37) 

For £11 = 1, the critical value of £12 is equal to 25. A comparison with the numerical solution of the GP equations 
shows that this value is overestimated: e.g., at M = 87, the efficiency T at the optimal time is ^ 0.9, while it should 
be equal to 1. However, the result ([37)1 gives a reasonable estimate of the point at which is no longer possible to 
transfer with perfect efficiency all the particles from one condensate to the other, due to the self-trapping of the initial 
condensate wavefunction. Finally, we observe that, for M < 0, the agreement between numerical and variational 
results is still good and the critical point is Mc = —47. We also notice that similar results can be obtained for the 
Rabi switch of iV-soliton solutions. 



B. Moving bright solitons and dark solitons 



The proposed protocol works also for transferring moving solitons: with £11 = £12 ~ £22 = one can prepare the 
initial wavefunction as 

^(^ 0) - c'- ( V^me'-^^'^ \ .33. 
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FIG. 3: Comparison of the population imbalance ri{t) obtained from the numerical solution of the GP equations (I25l) - (|26l) [solid 
line] with the results of the two-mode equations (|34|l [dotted line] for gi2 = —1,-13 given by Eq. (|36|l to a mass pendulum 
M = 0, 2, respectively. Inset: efficiency (solid line) and fidelity (dashed line), respectively defined according to Eqs. (|14[) 
and (|15p . and obtained from the numerical solution of the GP equations. Parameters used in both plots: to = 0, 7 = 1, 
gii = 522 = -1. 



For t <to one has (with fj. — —r/8) 



> " cosh {l{x - vt)/2) \ y]V2to)e^^^(o) J ' ^ ' 



so that at ti, i.e., at the end of the pulse, S = 7r/(27), one has 

.0 / / AT ^fr\\^iV 2(0) \ 

(40) 



coshii{x-vh)/2f ^ y]Vi(0)e^^i(o) 

If no particles are initially in |2), then one can transfer the moving soliton from a condensate to the other, as depicted 
in Fig. [Id)-(f). 

When the ^^j's are different, one can use the same variational method discussed in the previous subsection. One 
needs to consider the variational wavefunction (|3ip. but with a time-dependence included in the functions <i>, which now 
read ^i{x,t) = (l/2)e™^^™ cosh[£ii{x — vt)/2]: apart form constant terms, we obtain the same Lagrangian 

[c.f. Eq. (jB2p ] and the analysis is the same as before. In particular, the threshold for the self-trapping transition is 
still given by Eq. ([HT)) . 

On the other hand, when the parameters gij are positive and equal (and Vi = 0), the soliton solution is now a dark 
matter-wave soliton, and one can transfer it from one condensate to the other. To examine the situation when the gij 
are different and estimate the self-trapping threshold, one can also use a variational approach. Omitting the details, 
when gii = 322 one gets Mc — 47, with M = 311(312 — 3ii)'t-, where n is the (asymptotic) density of the dark soliton 
for very large x. 



C. Effect of the trapping potential 

First, we consider repulsive condensates, in the ideal situation where gn = gi2 = 322 = 1. An example of the 
realization of the Rabi switch for the ground state of the system is shown in the top panels of Fig. [D we have 
considered an harmonic trapping potential of the form V{x) = (1/2)0^3;^. In the numerical simulations of the GP 
equations, we use fl — 0.08, — 10 and 7 = tt/IO; hence, after t > 15, the condensate wavefunctions have completely 
switched between components. Next, we consider the attractive case with gn — gi2 — 322 = — 1- In this case, as an 
initial condition (for the first component) we have used a bright matter-wave soliton with the well-known sech-profile 
given by Eq. l|27|) . As shown in the middle panels of Fig. |4l the transfer of the wavefunction is complete, just as in 
the repulsive case. 

To better illustrate the versatility of the Rabi switch, even for A/" > 2, in Fig. [5] we have also considered the case 
with = 3. In accordance with the results of the analysis carried in Section III, we use ai = a2 = (7r/10)/V2 between 
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FIG. 4: Panel (a) shows a space-time plot of the density for the first component (top panel) and \ip2{x,t)f for the 

second component (bottom panel). Panel (b) shows the spatial profile of \Tpi{x,t — 0)|^ and |i/'2(a;,t — 25)|^ in solid and dashed 
lines respectively. The dash-dotted line shows the magnetic trap potential. They also show the evolution of the particle number 
Ni — f I'ipil^dx for each of the components in the interval of the dynamical evolution. These features are shown in panels (a) 
and (b) for gn — gi2 — 322 ~ 1. In panels (c) and (d), they are shown for a soliton in the case of gn — gi2 = 522 = — 1. 
Analogous features are shown for 3-component condensates in panels (e) and (f) (the third component is shown by dotted line 
in the right panel), again for the case with gij — 1 for all i,j = 1, 2, 3. 



to = 10 and t = 15. In line with Eq. ([23]) . at the end of this time interval, r2 = — 1 and, as a result (due to the 
symmetry in the choice of ai and 02), half of the matter initially at the first component is transferred to the second 
component and half to the third component, in excellent agreement with the results shown in the bottom panel of 

Fig. m 

When the gij's are different, the transfer will no longer be complete. As a measure of the deviation from the "ideal 
switch" , we characterize the degradation of the switch in this inhomogeneous case according to the relevant quantities 
introduced in Eq. and (|15p . For repulsive condensates, we have considered the case of the ground state of *^Rb, 





FIG. 5: Panel (a) shows the initial condition: all the particles are in the ground state of the first component. The solid line 
shows the harmonic trapping potential. The parameters used in the GP equations are fl = 0.08, to = 0, 7 = tt/IO, gu = 0.97, 
gi2 ~ 1.03. In (b) we plot the population imbalance r]{t) for gi2 ~ 1 (solid line), —1 (dashed line) and —10 (dotted line). The 
optimal times with the maximum transfer are respectively ti = 5,5.02,6.48. In (c), (d) and (e) we show the spatial profile 
p2{x,ti) = \4>2{x,ti)f at these optimal times. In (e), where the transfer is not optimal, we plot also the the spatial profile 
pi{x,ti) = \il)i{x,ti)\'^ . In panel (f) we show the transfer efficiency function T (solid line) and the fidelity F (dashed line) vs. 
the value of the inter-species interaction coefficient gi2 . 



where the two spin states mentioned above have corresponding scattering lengths such that gu : (?22 — 0.97 : 1.03. In 
this context, we have identified the ground state configuration for the first component alone and have subsequently 
apphed the linear coupling for 7 = tt/IO, for various values of 1712. The resulting transfer efficiency as a function 
of gi2 is shown in Fig. [5] (bottom right panel). The numerical result indicates that in the defocusing regime of 
repulsive inter-species interactions, the transfer efficiency and fidelity remain very high (> 0.9) throughout the interval 
— 10 < gi2 < 10. Similar results have been obtained for attractive interactions gu = 322 = — 1 and varying 1722, starting 
form the ground state of the first component in presence of the trap. 

We also performed similar computations in the case of attractive intra- and inter- species interactions for the 
solitonic initial condition (j27p in the first component, which, with the trap confinement, is no longer the ground state. 
We have set the intra-species interactions at 1711 = (722 = — 1, varying 1712. Our findings are plotted in Fig. [SI showing 
the robustness of our protocol in a range of values of 1712 between —3 and 3, a range smaller with respect to the case in 
which the initial condition is the ground state of the first component. This is due to the fact that the initial condition 
is not the ground state, and to the fact that the Rabi pulse used is homogeneous in space; note that (according to 
the analysis of Appendix A) a sort of space-dependent Rabi pulse would be needed to improve further the transfer 
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(a) (b) 




FIG. 6: Panel (a) shows the initial condition: all the particles are in the soliton (|27|l of the first component. The parameters 
used are Q = 0.08, to — 0, ^ = tt/IO, gu — gi2 = —1. In (b) we plot the population imbalance ri{t) for gi2 — —1 (solid line), 
1 (dashed line) and 3 (dotted line). The optimal times for maximum transfer are respectively ti — 5,5.10,5.58. In (c) and (d) 
we show the spatial profile of p2{x,ti) = \4>2{x,ti)f at these optimal times for gi2 = 1 and gi2 = 3. In panel (e) we show the 
transfer efficiency function T (solid line) and the fidelity F (dashed line) vs. the value of the inter-species interaction coefficient 
g\2, and in (f) we plot the same quantities for = —1, g22 ~ 1 and the same initial condition: the efficiency is reduced with 
respect to the case gu — g22 ~ —1. 



efficiency. Finally, as expected, the Rabi switch is less effective in the case in which gu and 522 have opposite signs. 
The quahtative reason is that the wavefunction, being a bright sohton which is supported by a negative coupling 
constant, cannot easily be transferred to an environment characterized by a positive coupling constant (which does 
not support bright solitons but rather dark ones). The fidelity and the efficiency are plotted in the panel (f) of Fig. 

El 

V. RESULTS FOR 2D SETTINGS 

Let us now consider the 2D version of Eqs. ([T])-([2]), pertaining to "pancake" -shaped condensates 0. We will focus 
on the realistic case of a binary mixture of two hyperfine states of ®^Rb, with gu : gi2 : 322 = 0.97 : 1 : 1.03, and 
examine both ground and excited states; the latter, will be characterized by the presence of one vortex or of many 
vortices arranged as vortex-lattice configurations. We shall consider only repulsive intra-species interactions, since 
for attractive ones, the system is generally subject to collapse |2|]. In order to compare the results with the ones 
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- y = n/lO 

- Y = Tt/20 
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§12 




FIG. 7: (a) The transmitivity function T versus the inter-species strength gi2 for 7 = n/10 (solid line) and 7 — n/20 (dashed 
line) . (b) Contour plot of the initial density of the ground state in the first component tpi . (c) Contour plot of the corresponding 
final density \4'2{t = 25) p of the second component for 312 — 1; here the transfer of matter is complete, (d), (e) Contour plots 
of the corresponding final densities of both species, namely \4>i{t — 25) ^ and \1p2it = 25) but for gi2 — 2; as it is seen part 
of the matter remains (is missing) in the first (from the second) component. 



pertaining to the ID setting, we vary the inter-species strength 1712 in the same domain (although for 1712 < the 
system is also subject to collapse). In a real experiment this may be done by using external magnetic fields, which 
can change the magnitude and sign of the scattering length through the Feshbach resonance mechanism f]\ . 

It is relevant to consider at first the efhciency T as a function of 1712, for different values of 7. The result is shown in 
the top panel of Fig. [71 In the 2D case, the transfer is almost complete as T > 0.95 for all positive values of gi2 < 2 
(for 7 = tt/IO). It is also seen that the efficiency is higher for larger values of the magnitude (or inverse duration) of 
the linear coupling coefficient 7. 

We have considered the ground state of the system, shown in Fig. El^b), in the first component ^i, for an harmonic 
trapping potential with strength Q = 0.045 (the chemical potential is equal to one). In this case, assuming that the 
switch parameters are to = 10 and 7 = tt/IO, we have found the following: for gi2 = 1, the transfer of matter in 
the second component ?/'2 is complete [see Fig. [7{c) where the density of the second component at i = 25 is shown], 
while for 312 = 2 it is incomplete. In particular, as shown in Figs. [7|^d)-(e) (where the densities \ipi{t — 25)p and 
\4'2{t — 25) p are respectively shown), a fraction of matter remains in the first component and is correspondingly 
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FIG. 8: (a) Same as in Fig. [T^a). (b) The initial density of the first component ipi, consisting of a cloud with one vortex in the 
center, (c) The corresponding final density \'ip2{t = 25) |^ of the second component for gi2 — 1; here the transfer is complete 
and the final configuration is identical to the initial one. (d), (e) The corresponding final densities of both species, namely 
\ipi{t = 25)p and \^2{t = 25)^, but for gi2 = 2. 



missing from the central part of the second component after the switch-off of the Rabi pulse. 

Next we consider an excited state, in which a vortex is initially placed at the center of the BEC cloud (first 
component). Here, it is interesting to investigate whether such a coherent nonlinear state can be transfered in the 
second component. As seen in Fig. [5l^a), the efficiency is as high as for the ground state transfer. Also, for gi2 = 1, a 
perfect transfer of this excited state occurs as seen in Figs. [SJb) and (c), where the initial state of the first species and 
the final one of the second species are respectively shown. Nevertheless, for gi2 — 2, the transfer is not complete, as 
it can be seen in Figs. EUd) and (e): starting again from the initial density shown in in Fig. Efb), after the switch-off 
of the process, a ring-shaped part of the matter remains in (is missing from) the first (second) component. A careful 
observation of Fig. [SJd) also shows that this "high" density ring surrounds a low density (« 0.2) part of matter with 
a vortex on top of it. Note that even in this case, the vortex is transferred in the second component; on the other 
hand, the above mentioned "bright" and "dark" ring structures (respectively in the first and second component) do 
not carry any topological charge. It is interesting to remark that such methods are similar in spirit to the ones used 
to produce ring-like patterns in the recent experimental work of [20j . 

Finally, we have considered a vortex cluster, namely a triangular vortex lattice, initially placed on top of the BEC of 
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FIG. 9: (a) Same as in Figs. [7ja) and[8ja). (b) The initial density of the first component t/ji, consisting of a TF-cloud with 
a triangular vortex lattice (« 24 vortices), (c) The corresponding final density \Tp2{t = 25) |^ of the second component for 
gi2 = 1; here the TF cloud and the vortex lattice are perfectly transferred to the second species, (d), (e) The corresponding 
final densities of both species, namely \ipi{t = 25) ^ and \ip2{t — 25) p, but for gi2 — 2. 



the first component. Similarly to what was found for the ground state and the single vortex, we find that the transfer 
efficiency function shown in Fig. [DJa) assumes values very close to 1 for a wide range of values of gi2 (especially so 
in the case of short pulse durations i.e., fast transfer). As shown in the example of Figs. [Hb) and (c) for 512 = 1 the 
vortex lattice in the first component is perfectly transferred to the second one. An "imperfect transfer", for gi2 = 2, 
is shown in Figs. O^d) and (e), corresponding to the final states of the first and second species when the initial state 
is as in Figs. [DJb). Here it is observed that small spots of matter (with densities « 0.025), each of which carries a 
vortex on it, remain in the first species, while a robust crystal structure is transferred to the second component. It is 
worth noting here that the vortex lattice remarkably appears to be even more robust, in its "switching properties", 
than the ground state of the 2D system. 
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VI. CONCLUSIONS AND OUTLOOK 



In this paper, we have systematicaUy examined the possibihty of using the hnear Rabi couphng between the two 
(or more) components of a Bose-Einstein condensate as a means of controUably transferring the wavefunction of 
one condensate to the other. In particular, we have focused on the case of different hyperfine states, even though 
similar considerations are applicable to multicomponent condensates composed by different atomic species. We have 
illustrated that this transfer is exact and can be analytically studied in the limit where all inter- and intra- species 
interactions are equal. In addition, we have studied departures from this limit both numerically and by means of 
a two-mode ansatz, showing that in this two-mode description the impossibility to transfer all the particles from a 
condensate to the other is seen as the self-trapping of the initial condensate wavefunction. The threshold for the self- 
trapping has been compared in the homogeneous limit with the findings of the numerical simulations Gross-Pitaevskii 
equations. 

The two-mode analysis shows that that for deviations from the ideal case one can optimize the transfer by choosing 
an appropriate pulse duration different from the duration given by Eq. (jlip for equal interaction strengths . In the 
presence of external trapping potentials, our numerical simulations show that, for repulsive condensates (but, to a 
lesser extent, also for attractive condensates), the Rabi switch is very robust with high efficiency. In general, a change 
of sign of the inter-species interactions in ID has been shown to degrade -although, under appropriate conditions, not 
substantially- the efficiency of the transfer process. The switching was surprisingly found to be even more robust for 
particular types of coherent patterns, such as vortex lattices in pancake-shaped condensates. Furthermore, we have 
also illustrated that the generalization of the proposed Rabi switch to more than 2 components can offer a possibility 
for systematically routing matter, in a controllable way, between different "atomic channels" . 

It would be quite challenging to examine experimental realizations of such atomic switches and routers especially 
within the context of two hyperfine states, but also in that of spin-1 states recently studied in [l^, [l3]- Of particular 
interest in this setting would be the dynamics and phase separation of vortices and vortex lattices in multicomponent 
condensates. We expect that our mechanism may be relevant when transferring a condensate wavefunction also in 
different components of a dipolar BEG [60| . 

Another relevant issue concerns the analysis of the role of the quantum fluctuations on the Rabi switch of wave- 
functions: indeed, for two condensates (whose dynamics is described by coupled Gross-Pitaevskii equations) the Rabi 
switch is not copying the full many-body wavefunction of a component in the other, but only the (one-body) or- 
der parameter. Then, we expect that the presence of quantum fluctuations would naturally degrade the efficiency 
of our protocol. A study of such degradation is relevant for the implementation of a quantum register through a 
two-component Bose gas in an array of double- wells. 

The protocol proposed in this paper allows the possibility to copy the wavefunction of a condensate into another 
species, providing a matter- wave counterpart for optical switches realized in nonlinear fiber optics [l^]. "All-optical 
switches" are very important tools in this context, allowing for to manipulate and eventually store the information 
contained in optical solitons. For this reason, the study of protocols focused on copying of wavefunctions in Bose 
condensates is deemed to be an important task towards an efficient manipulation of matter-wave solitons. Furthermore, 
one can think of a possible application of the present protocol to implement copies for quantum registers: indeed, one 
could perform some operations on a component (e.g., with a component in an array of double wells) and arrive at a 
desired wavefunction. At this point the wavefunction can be copied on the other component, and one can gain access 
(at a later time) to this information. In this respect, it would be very important to investigate the degradation of the 
quantum switch due to the quantum fluctuations. 

Finally, since the protocol proposed in this paper is not restricted to linearly coupled Gross-Pitaevskii equations, 
one may safely expect that - once the effects of quantum corrections are accounted for - it can be extended to other 
sets of coupled equations such as the ones arising in the analysis of the weak pairing phase of p-wave superfluid states 
of cold atoms [6l| and of quantum wires embedded in p-wave superconductors [65| . 
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When the interaction strengths gij are different (and, for simplicity, Vi = V2), then Eqs. ©-(HI) can be written as 



APPENDIX A: DIFFERENT POTENTIALS AND INTERACTION STRENGTHS 




(Al) 
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with 

L)^ <^^' 

Upon removing the Rabi term in Eq. (jAip by setting -0 = U(f), with U given by Eq. one gets 

»^ = -^AV; + V0 + ^(0^Gj0K-0, (A3) 
i=i 

where Gi = £i — 5^(£i — £2) and G2 = £2+'5^(£i — £2), where the 2x2 matrices £1^2 are defined by £1 = Gi — iSSiA 
and £2 = G2+iSS2A, with A ~ —iS(a\~ 02) -^C{aj, — a^, S = sinX(<), C = cosl(t) and Si = 511 — 512, <52 = 522 —512- 
Furthermore G3 = —iCS{Ci — £2) = — G4. For 511 = 512 = 522 = 5, then 61 = 82 = Q and Gi — G2 — G and 
G3 = G4 = 0, so that Eq. ([7]) is retrieved. 

With different external potentials {Vi{r) ^ V2{r)) and equal interaction strengths (511 = 522 = 512 — 5), Eqs. 
([T])-© can be written in a matrix form as 

^^ = \ (-iv)% + (i^tG^) ^ + Vi{Y)ai^ + V2{Y)a2^ + a{t)Pi,, (A4) 

where 

= ( ) ' = ( E! ? ) , ^3 - ( J ) - ( J ) , (A5) 

and G given in Eq. (|4]). It is yet possible to perform a decomposition permitting to formally write Eq. (|A4p without 
the Rabi term a{t)P'ip. We set ip — U(f>, with 

C/ = f "1 "0 : (A6) 

\Ui U2 J 

the Rabi term vanishes, provided that the functions Uj{r,t) (j — I, - ■ ■ ,4) obey the (linear) matrix equation 

= -^AU + a{t)PU + V{r) {u^a^ - mcj^) , (A7) 

where V{r) — Vi{r) — V2{r). This matrix equation corresponds to four equations for ui, U2, "a, and U4, which are 
grouped in two pairs, one for ui and M4, namely 

= -2^1*1 +a(t)M4, (A8) 

«^ = --Au4-V(f)u4 + a(t)wi, (A9) 

and a similar for U2 and U3 (with V(r) instead of — V(r)). Eqs. (|A8|) - (|A9p are two coupled linear Schrodinger equations, 
and the difference between the potentials V enters as an effective potential in one of the two equations. When V = 0, 
the result ui — U2 — cosX(t) and 1/3 = ^4 = — isinX(i) is readily obtained. 
With the functions Uj defined by Eq. (|A7[) . then (j) satisfies the equation 



i (-zv)^0+ ((/.tG^) <?!) + Vi(r)cri0 + y2(r)a2(/> + A'- V(/), (AlO) 
with G = WGU and X = X{r,t) ^ -U-^WU. Eq. ((IlO| can be written in the form 

i (-iV + i-Y)%+ (0^G0) <?!. + yi(r)ai0 + V2(r)cr20 + y(/), (All) 



where y = yX^ — V • Xj /2, showing that, with different external potentials, an effective complex vector potential 
acts on the multicomponent gas and the effective interaction strengths are both time- and space- dependent [since 
the G depends upon the U{f, t)]. 
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APPENDIX B: TWO-MODE VARIATIONAL EQUATIONS 



For the pulse ([8]) and at times to < t < to + S, the Lagrangian to be computed is given by £ = ^{ipy^^^ — —g^i^v) — 
{ipyHtpv), where 

_____ ^ o2 

V'y is defined in Eq. (|3T|) and iiT = ~2'§^- Then, one gets 

C = -m^i - N2^2 - l\h^F f ^) ^Njhcos (ipi - ^2) - - %-N2 + ^l+ 

V til V^ii/ ^4 z:4 

where the functions F{d) and £(6*) are defined as 

Fie) = r Lie) = r 

cosh (y) cosh (%)' y_oo cosh^ (y) cosh' (%) ■ 

Note that F(0) = tt, F(1) = 2, L(0) = 2, and L(l) = 4/3, while both functions ^ as ^ 00. For smaU values of 

56 ^e -I, one has F(6l) w 2 - (561 + "^^^^e"^ and L(6') w | - fje* - ^^5^56*2. For £11 = £22, from Eq. (jB2| one 
gets, apart from constant terms, the Lagrangian 

C = -iVi<^i - iV2(^2 - 27y]ViiV2 COS (^1 - if 2) + ^A^f + IfA^I + ^-^4^N^N2. (B4) 

12 12 D 

Introducing the variables p3p . the equations of motions for r/ and ip obtained from the Lagrangian (jB2p are 



?7 = 27' ^1 — rf sin(y9, 

(yj = -2i^ — cos ip + /\E + Mr], (B^) 

^l-J)2 



where 

7' = ,, 

2 V ^11 \«ii 



7 ("^Ei ) , (B6) 



M = — 
12 



2 V^ii 



^2 _ /i2 
■'^11 *22 



and 



1 



The variables ry, are canonically conjugate dynamical ones with respect to the effective Haniiltonian 

Heff = + 27V1 - "if cos(y5 + A£; • 77. 



(B7) 



(B8) 



(B9) 



Equation (|B9p is the Hamiltonian of a non-rigid pendulum, whose mass and length depend on the parameters iij. 
For £11 — £22, then M = £11(^12 — ^ii)/6, AE = and 7' = 7, retrieving Eq. (f34|) and the effective Hamiltonian 
([35]) . When £12 = ^11 — £22, then the mass of the pendulum is vanishing. When in ^ £22, then the presence of the 
detuning AE favours self-trapping and can be studied as in (STj . 
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